Cosmic string-seeded structure formation 
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We describe the results of high-resolution numerical simulations of string-induced structure for- 
mation in open universes and those with a non-zero cosmological constant. For models with 
r = Qh = 0.1-0.2 and a cold dark matter background, we show that the linear density fluctuation 
power spectrum has both an amplitude at 8/i~^Mpc, as, and an overall shape which are consistent 
within uncertainties with those currently inferred from galaxy surveys. The cosmic string scenario 
with hot dark matter requires a strongly scale-dependent bias in order to agree with observations. 
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A. Introduction — In this Letter we describe new re- 
sults from an investigation of cosmic string-seeded struc- 
ture formation in hot and cold dark matter models. The 
cosmic string scenario W| predated inflation as a realistic 
structure formation model, but it has proved computa- 
tionally much more challenging to make robust predic- 
tions with which to confront observation. The present 
paper relies on high resolution numerical simulations of 
a cosmic string network M with a dynamic range extend- 
ing from before the matter-radiation transition through 
to deep in the matter era (developing on previous work 
[g[). We calculate the linear power spectrum of density 
perturbations V{k) induced by the strings in flat models 
with and without a cosmological constant, and we then 
extrapolate to open cosmologies. This work represents a 
considerable quantitative advance by incorporating im- 
portant aspects of the relevant physics not included in 
previous treatments. 

In the first instance, we consider density perturba- 
tions about a fiat Friedmann-Robertson- Walker (FRW) 
model with a cosmological constant A and which are 
causally sourced by an evolving string network with 
energy-momentum tensor 0c(^(x, 77). In the synchronous 
gauge, the linear evolution equations of the radiation and 
cold dark matter (CDM) perturbations, d^ and S^ respec- 
tively, are given by (modified from M) 
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where 84. = 800 + 8^^, a is the scale factor, the sub- 
script "eq" denotes the epoch of radiation-matter den- 
sity equality, "0" denotes the epoch today, a dot rep- 
resents a derivative with respect to conformal time 77, 
r^c — 8ttGpco/3Hq and ft a — A/3Hq. It proves use- 
ful to split these linear perturbations into initial (I) and 
subsequent (S) parts [||, 5n{x, 77) = Slf{x, rf) -f 5%{'x., rf) , 
where N = c, r. The initial perturbations (5^(x, 77) de- 
pend on the string configuration at some early time rji, 
because the formation of strings creates underdensities in 
the initially homogeneous background out of which they 



are carved. The subsequent perturbations (5^(x, 7;) are 
those which are generated actively by the strings them- 
selves for 77 > 77i. Because strings induce isocurvature 
perturbations, (5^(x, 77) must compensate (5^(x, 77) on co- 
moving scales |x — x'l > 77 to prevent acausal fluctuation 
growth on super horizon scales. 

The system of equations (|l|,0) can be solved for the 
subsequent perturbations 5^ (x, 77) , with initial conditions 
(5f = (5f = and jf = jf = at 7/ = ryi , by using a discretized 
version of the integral equation with Green's functions: 

4(x,77)=47rG/d7/ /dVgAr(A;7/,r/)8+(x',77'), (3) 



where A = |x — x'|. The Green's functions in Fourier 
space can be calculated numerically by solving the homo- 
geneous version of (|]J|) with initial conditions at 77 = 7;': 

^c = 3^,/4=landge=^r = (^Ar = for 77<77'). 

The subsequent perturbations 5'^ (x, 77) are dynamically 
sourced by moving local strings with spacetime trajecto- 
ries we can represent as a;^ — (77,Xs(i7, 77)), where tr is a 
spacelike parameter labelling points along the string (a 
prime represents a derivative with respect to a). The 
stress energy tensor of the string source is then given by 
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where /i is the string linear energy density, e = [x^ /(I — 
Xg)]^/^, and we have also assumed that Xg • Xg' — 0. In 
this case, it is straightforward to compute 8+ in (Q) as 



8+(x,7;)=8 
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The stress energy 8^,y was calculated directly from high 
resolution string network simulations [p[. Dynamical 
ranges exceeding 100 in conformal time (redshifts up to 
1000) were achievable because of a 'point-joining' algo- 
rithm maintaining fixed comoving resolution [p| and par- 
allelization. 

B. Approximation schemes — It is a very substantial 
numerical challenge to evolve the initial and subsequent 



perturbations induced by cosmic strings such that they 
accurately cancel on superhorizon scales by the present 
day rjQ. For the large dynamic range required for the 
present study, we have by necessity adopted the 'compen- 
sation factor approximation' suggested in a semi-analytic 
context in ref. H. To implement this, we accurately 
evolved the long string network numerically — the domi- 
nant active source term — and then multiplied the Fourier 
transform of the resulting stress energy 9+(fc,77) by a 
cut-off function F{k,7j) = [1 + (fcc/A;)^]"^ This results 
in the correct k'^ fall-off in the power spectrum at large 
wavelengths above the compensation scale k~^ r^ rj. The 
efficacy of this approximation has been demonstrated by 
studying multifluid compensation backreaction effects in 
ref. 0. For the present study we have adopted the an- 
alytic fit for kc{ri) presented in ref. 0, which smoothly 
interpolates from kc = \Qrj~^ in the radiation era to 
fcc = vlSr;^^ in the matter era. The quantitative im- 
plementation of compensation is a subtle issue and a key 
uncertainty in all work to date on gauged cosmic strings. 
We note, however, that our results are relatively insensi- 
tive to the choice of fee, especially in open and A- models 
(e.g. a large factor of 2 increase in k^ causes only about a 
20% decrease in the power spectrum at fc « 0.15/i~^Mpc 
for a flat J^a — 0.8 model). 

In order to study the formation of structures with cos- 
mic strings in hot dark matter (HDM) models, we use a 
reasonably accurate alternative to much more elaborate 
calculations using the coUisionless Boltzmann equation. 
We simply multiply the Fourier transform of the string 
source term 0_|_(fc,?7) by a damping factor G{k,ri) — 

[1 + (0.435fcL»(?7))^"^]~'*'*^ g. Here, D{t]) is the comov- 
ing damping length that a neutrino with velocity T^/rrii, 
can travel from the time rj onwards. The factor G{k, rj) is 
a fit to numerical calculations of the transfer function of 
a Fcrmi-Dirac distribution of non-rclativistic neutrinos 
and accounts for the damping of small-scale perturba- 
tions due to neutrino free-streaming [p|. We calculated 
D{ri) numerically and found an excellent fit to our results 
forT^o = imUxlQ-^^GeY andm^ = 91.5 rj/i^eV with 

The other key difficulty confronting defect simulations 
is their limited dynamic range. At any one time, an 
evolving string network sources significant power over 
a lengthscale range which exceeds an order of magni- 
tude. However, we can employ a semi-analytic model 
to compensate for this missing power |q|, which proves 
to be fairly accurate in the scaling regimes away from 
the matter-radiation transition. The procedure is es- 
sentially to square the expression (y) in Fourier space 
to obtain the power spectrum V{k). This becomes a 
Green's function integral over the unequal time corre- 
lators (0+(k, ?7)8+(— k, 77')), which are subsumed in a 
scale-invariant 'string structure function' JF(k, 77), that 
is, the power spectrum is given by |o]: 



V{k) =. levr^GV' / \gN{k;m,ri)\''T{k,ri)din. 

Jrji 
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In practice, we obtained the structure function J^(k, 77) 
phenomenologically by fitting its shape and amplitude to 
simulations of limited dynamic range deep in the matter 
and radiation eras. We were then able to use an interpo- 
lation based on the actual string density during the tran- 
sition era to provide a good fit to the simulation power 
spectrum for any given dynamic range ryi — > 77. 

C. Open and A-cosmologies — The possibility that the 
universe is open or has a cosmological constant is now 
favoured by a number of oberservations, so it is natural to 
explore the string-induced spectrum in these two regimes. 
For the A- models ilc + I^a = 1, we have taken the weak 
dependence for the COBE normalisation of Gn ex 51""°^ 
suggested in ref. ||] for (||,||). For the open models fie < 1 
and A = 0, we simply rescale the simulated spectrum 
from a flat universe with A = in the following way 
(adapted from ||]): 

s{k, h, n,) = s{k, 1, 1) • nlh^ ■ f{n,) ■ g^{n,), (7) 

where k is in units of Qch^ Mpc^^ and f{flc) = ^^cT"'^ I'fi- 
flects the COBE normalisation of Gfj, ||]. The last factor, 

g{^c) = 2.5ric/(l + ^c/2 -I- ric ), gives the total sup- 
pression of linear growth for density perturbations in an 
open universe relative to an flc~l and J^a = universe 
prt. Finally, the second factor fl^h'^ in (H) arises natu- 
rally from the normalization of the Green's functions in 
(0,0). We have verified that a similar analytic rescaling 
from an r^c — 1, ^A = model to Qc + ^A = 1 rnod- 
els agrees very accurately with the solutions of (|l|,0) ob- 
tained from simulations. This also helps to justify the 
extrapolation still required for open models, although 
some uncertainty remains concerning the COBE normal- 
ization. 

D. Results and discussion — In figure |l| we plot the 
Cosmic Microwave Background (CMB) normalized (i.e. 
GfiQ = Gfi X 10® = 1.7 [11|) linear power spectrum in- 
duced by cosmic strings in an Oc = 1 CDM cosmol- 
ogy with h = 0.7. The central set of numerical points 
was sourced by string network simulations beginning at 
rj = 0.47yoq which were continued for 1318 expansion times 
(from redshift Zi ~ 31700 to Z[ « 23). String simula- 
tions were always ended before the horizon grew to half 
the simulation box-size. These had a string sampling 
resolution at least four times higher than the structure 
formation grids, which had up to 256"^ points with over- 
all physical scales ranging from 4-100ft-~^Mpc. Given the 
dynamic range limitations, we have also plotted the semi- 
analytic fit m) over the full range of wavcnumbcrs, illus- 
trating the good agreement with our numerical results. 
This was also apparent in HDM simulations, so we have 
considerable confidence that the semi-analytic model pro- 
vides a good approximation to the shape and amplitude 
of the string simulation power spectrum. These results 
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FIG. 1. Comparison of the string-induced simulation power 
spectrum and the semi-analytic fit (n). The top-right, central, 
and bottom-left solid lines with crosses are the simulation re- 
sults in the deep radiation, transition, and deep matter eras 
respectively. The dashed lines are the semi-analytic fits cor- 
responding to the same dynamical ranges of the simulations. 
The solid line is the semi-analytic model over the full dynamic 
range from ?7i = to today. 



are also qualitatively consistent with the semi-analytic 
results of re f. [p[ and also with unpublished matter era 
simulations ]12( |. 

In figure H, we make a comparison between our CDM 
and HDM string power spectra and the observational re- 
sults inferred from galaxy surveys |l3], in each of five 
different background cosmologies: (I) 51c. h = 1, f^A = 0, 
(II) flc,h ^ 0.3, f^A = 0.7, (III) f7c,h = 0.3, 17a = 0, (IV) 
flc,h = 0.15, VLk = 0.85, and (V) fic,h = 0.15, Oa = 0. 

Consider first the ilc = 1 CDM model. We calculated 
the standard deviation of density perturbations a% by 
convolving with a spherical window of radius 8h~^ Mpc 
to find cr8(sim)(^ = 0.5) = 0.32Gne, (Ts(sim)ih = 0.7) = 
O.SQG^e and as(^si^^{h — 1.0) — OAlGfiQ. A compari- 
son with the observational data points shows that strings 
appear to induce an excess of small-scale power and a 
shortage of large-scale power, that is, the Oc = 1 string 
model requires a significant scale-dependent bias. This 
is not necessarily a fatal flaw on small scales because, as 
the corresponding HDM spectrum indicates, such excess 
power can be readily eliminated in a mixed dark matter 
model. However, the problem is less tractable on large 
scales where biases up to crioo(obs)/cioo(sim) ~ 3.9 around 
100/i~^Mpc might be inferred from the data points (using 
G/i6 = 1.7 and h = 0.7). Should we, therefore, rule out 
string models on this basis fljl? Although the flc = 1 
spectrum looks unattractive, there are three important 
mitigating factors. First, the present observational de- 
termination of the power spectrum around 10Gft,~^Mpc 
is very uncertain. Secondly, the immediate nonlinearity 
of string wakes means that strong biasing mechanisms 
might operate on large scales. Finally, unlike inflation. 
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FIG. 2. The CMB normalized linear spectra induced by 
cosmic strings in CDM (left) and HDM (right) models for 
different background cosmologies. The solid lines are the fiat 
A-models; the dashed lines are the open models. Here, we 
use h — 0.7 and the data points with error bars are the linear 
spectrum reconstructed from observations [[13[. 
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FIG. 3. A comparison of the observationally inferred mass 

fluctuation crgjobs) with that induced by cosmic strings in our 

simulations, o-gjsim). o-8(obs) is shown as the shaded area [jl5|, 

while a8(sii„) is plotted as dot-dashed {h = 1.0), solid {h = 0.7) 

and dashed(/i = 0.5) lines. 



defect models have never been wedded to an fl = 1 cos- 
mology. 

[b] We can observe from figure g, that for open or 
A-cosmologies with ilc « 0.1-0.3, the string + CDM 
power spectrum is much more encouraging. We find 
that the bias on large scales is always reasonably close 
to unity and, overall, it is much less scale-dependent. 
For example, over the full range of lengthscales in model 
IV (fie = 0.15, Oa = 0.85), the relative bias remains 
CTioo(obs)/c^ioo(sim) ~ 1.4 ±0.2 at 100/i-^Mpc. In figure 
S, the value of <78(si,„) induced in our simulations with 
the CMB normalized Giiq P,|ll[ is compared with the 
observationally inferred cr8(obs) |15| for the full gamut 
of open and A-models. We can see from figure ^ that 
f^8(obs)M(sim)(f^c=l)« 0.79 ±0.21, 0.95 ±0.25, 1.17 ± 
0.31 for h = 1.0, 0.7, 0.5 respectively. When h = 0.7, 



o'8(sim) matches cr8(obs) within the uncertainties for flat 
A-models when J7c > 0.35 and for open models when 
ilc > 0.4, while for both cases the ratio cr8(obs) /c8(sim) ^ 
2 for all flc > 0.1. Combining these results with an anal- 
ysis similar to figure 0, we found that the best string 
models lie in the range T = ilh = 0.1-0.2, producing 
both an acceptable cr8(sii„) and power spectrum shape. 
Hence, an open or A-cosmology in the context of string 
+ CDM model certainly merits a more detailed nonlin- 
ear study. These conclusions are in qualitative agreement 
with semi-analytic results P,|l6ll and those based on a 



phenomenological string model IH 

As for the HDM results, the comparison with observa- 
tion seems to require a strongly scale-dependent bias for 
any choice of the cosmological parameters (models I-V) . 
However, the lack of small scale power may be partially 
overcome if baryons are properly included in the analy- 
sis. Further investigation using a hydrodynamical code 
will be required to determine whether galaxies can form 
early enough. 

A key feature of all these string-induced power spec- 
tra is the influence of the slow relaxation to the mat- 
ter era string density from the much higher radiation 
string density, which has an effective structure function 
^{k, rj) in ra) with approximately 2.5 times more power 
than the matter era version. Even by recombination in 
an r^c = 1 cosmology, the string density is more than 
twice its asymptotic matter era value to which we nor- 
malize on COBE scales. This implies that the string 
model provides higher than expected large-scale power 
around 100ft-~^Mpc and below. Interestingly, this can 
also be expected to produce a signiflcant Doppler-like 
peak on small angle CMB scales, an effect noted in 
ref. ]ll| but not observed because only matter era strings 
were employed. Recent work in ref. [|l7| conflrms that 
such Doppler-like features can result from signiflcant non- 
scaling effects during the transition era. 

Finally, we comment on the fact that the key uncer- 
tainties affecting these calculations primarily influence 
the amplitude of the string power spectrum, rather than 
its overall shape which appears to be a more robust fea- 
ture. These uncertainties mainly result from the compen- 
sation approximation (mentioned previously) , the COBE 
normalization of the string energy density |0,nq| , the an- 
alytic approximation to the Green's functions, and sys- 
tematic errors [Q. Combining our best estimates of 
these uncertainties gives an approximate factor of 2 un- 
certainty in the power spectrum amplitude for f^c = 1 
and A-models. The extrapolations required for open 
models with F « 0.15 increase this uncertainty to at 
least a factor of 3 overall, but we will discuss this at 
length elsewhere [091. 

E. Conclusion — We have described the results of high- 
resolution numerical simulations of structure formation 
seeded by a cosmic string network with a large dynami- 
cal range taking into account, for the flrst time, modiflca- 



tions due to the radiation-matter transition. Our results 
show that for F = flh = 0.1-0.2 both as and the power 
spectrum shape of cosmic string-induced CDM fluctua- 
tions agree satisfactorily with observations. In particular, 
the generalization to open or A-models tends to remove 
the excess small-scale power found in cosmic string mod- 
els with ilc = I and SIa = 0, while also bolstering the 
large-scale power. The HDM power spectrum requires 
a strongly scale-dependent bias either on small or large 
scales, but we note that a high baryon fraction may help 
to increase small-scale power. We conclude that the pic- 
ture which emerges for particular cosmic string models 
seems encouraging and certainly deserves further study 
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